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ABSTRACT 

We present a deep, wide-field photometric survey of the newly-discovered Hercules dwarf spheroidal 
galaxy, based on data from the Large Binocular Telescope. Images in B, V and r were obtained with 
the Large Binocular Camera covering a 23' x 23' field of view to a magnitude of ^25.5 (5a). This 
permitted the construction of colour- magnitude diagrams that reach approximately 1.5 magnitudes below 
the Hercules main sequence turnoff. Three-filter photometry allowed us to preferentially select probable 
Hercules member stars, and examine the structure of this system at a previously unattained level. We 
find that the Hercules dwarf is highly elongated (3 : 1), considerably more so than any other dSph 
' satellite of the Milky Way except the disrupting Sagittarius dwarf. While we cannot rule out that the 

unusual structure is intrinsic to Hercules as an equilibrium system, our results suggest tidal disruption 
as a likely cause of this highly elliptical structure. Given the relatively large Galactocentric distance of 

(— I ' this system (132 ± 12 kpc), signs of tidal disruption would require the Hercules dwarf to be on a highly 

O |. eccentric orbit around the Milky Way. 

q , Subject headings: galaxies: individual (Hercules dSph) — galaxies: kinematics and dynamics 

&\ 

1. introduction In this paper, we present and analyse deep images of the 

i i ' a, ,i i , i • i p n i f i recently discovered Hercules dSph, which lies at a distance 
At the lowest luminosity end ot the realm ol galax- £ ^ i - — ni^nyh ggiok + & \ • 

ies, dwarf spheroidal (dSph) galaxies are characterised . ~, ,. p „" — 77 — ur . V ir , , — , U ^ V , , ' m ~ 

7— I ■ , ^ -1 • -i , x l-i ierred that the Hercules dSph displays an extended mor- 

>, as low surface brightness systems which are highly dom- . . . , . ,,. , , , ,. A 

,,,,, , , r pnology and may contain multiple stellar populations. A 

• mated by dark matter and have star formation his- , • , r77- T „ , — 1 , 

v-A , . c_ ,1 • 1 , , spectroscopic survey by Simon & Gcha 2007) measured a 

VO tones ranging from the simple to the complex. fhe ^ , ,• , , , ■ I ■ ; — 1 r 1 

^£ A/r-n ttt n.mn .1 v. • ^ 1 jc 1, mass-to-hght ratio ol 332±221 m solar units, and found ev- 

u") Milky Way (MW) currently has eighteen known dSph .. _ ° . ,. . T . ' 

■ c • 1 j- + i, tt 1 idencc lor kinematic substructure. In this letter, we exam- 

— 1 companions, many ol these, including the Hercules . ,. , , . ,. TT . irl . . ' 

, r , . , , , i-r, , \ — — 1 — f lonnTh ii me the structure 01 the Hercules dSph using photometry 

dwarf spheroidal galaxy Bclokurov ct al. 2007) recently ... . ... . , / . . br n , , J 



j , ,v o, tv ., , TT^ c /nncc obtained with the newly commissioned Large Binocular 

O discovered by the Sloan Digital Sky Survey (SDSS: „. TTr , ,,,,,,, , , r TT 1 

. nm at r< — Tu 1 — Tl[ori7v7\ ci i r 1 1 1-1 lelescope. We demonstrate that the structure of Hercules 

i~ — , Adelman-McCarthy ct al. 2007). Each ol these must be , , , „ , 



1 j ■ j • ' j ,1 ,„„ is highly elongated, much more so than all other known 

experiencing some structural distortion due to the MWs ,. . ,„ , . , . , , ,. , 

., , . i c i 1 i , • j 1 , ii- , ■ distant dSphs, possibly indicating strong tidal distortion, 

gravitational field, and tidal distortions and disruption fit- j t> t> 

have been recognized in some of them: foremost the Sagit- 

tarius dSph dlbata. Gilmore. fc Irwinlfl^l) but also the 2 ' DEEP photometry from the let 

Ursa Minor (|Martme~ Dclgad o et al.l [2001k iPalma et al.l The Large Binocular Telescope (LBT) is located on 



o 

> 

X 

I 120031) dSph. However, tidal distortion has not yet been Mount Graham in Arizona, a nd consists of tw o 8.4 metre 
observed in a satellite > 100 kpc from the Galactic centre, mirrors on a common mount (Hill et all 12006ft . Our data 
where the tidal forces should be much weaker. However, were obtained as part of the LBT Science Demonstration 
it is the pericentric distance of a satellite, not its current Time during which a single mirror of the LBT was fit- 
distance, that determines how much a satellite is affected ted wi th the blue channel of the Large Binocular C amera 
by tides. (XBC; iRagazzoni et al.l [20061 iGiallongo et al.ll2007D . The 
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LBC is a wide-field imager which provides a 23' x 23' field 
of view, sampled at 0.23 arcsec/pixel over four chips of 
2048 x 4608 pixels. LBC-Blue is optimized for the UV- 
blue wavelengths, from 320 to 500 nm, and is equipped 
with the U, B, V, g and r filters. 

The Hercules dSph was imaged with the LBT/LBC-Bluc 
setup on March 17 and May 13 2007. The data consist of 
five 5 min exposures in Gunn r-band, four 5 min expo- 
sures in l/-band and six 5 min exposures in B-band (all 
dithered), taken during photometric conditions with a see- 
ing of approximately 0.8 — 1.1". The data were reduced us- 
ing standard IRAF routines in the mscred package: the im- 
ages on the four chips were trimmed and bias-subtracted, 
and they were then flat-fielded using the combined twi- 
light flats obtained at the start of the night. The images 
contained significant field distortion at the edge, resulting 
from the 'fast' focal ratio (//1.14) of the LBTs primary 
mirrors. This was removed before co-adding the science 
frames, by applying a quadratic radial correction (accu- 
racy 0.2"). The Hercules images were then media n com- 
bined using the mscstack and mscskysub routines (|Valdesl 
2002]) to produce the final science frames. 

Photometry was measured using the PSF-fitting algo- 
rithm within the DAOPHOT package. The photometric 
dataset was calibrated using stars from SDSS, bypassing 
any need for an atmospheric extinction correction. The 
zero-point uncertainties were Sr ~ 0.02 mag and SB, 
SV ~ 0.03 mag, in cluding the empiric al colour transforma- 
tions described bv lJordi et alj (|2006l ). without any indica- 
tion of zeropoint gradients across the LBC field. We then 
correcte d the photometry fo r extinction using the dust 
maps of lSchleeel etU] (| 1998ft . with (E(B - V)) of 0.055 
mag. The photometric uncertainty and completeness as 
functions of magnitude were assessed using artificial star 
tests. We placed 1600 artificial stars in the image and 
attempted to recover them with DAOPHOT, where the 
photometric uncertainty was then determined as the dis- 
persion of the returned magnitudes about the mean (that 
is, not the input). This was repeated for artificial stars at 
every 0.25 magnitudes in the B, V and r frames. 

The resulting CMDs for selected regions of the image 
are shown in Fig. [1] For our analysis of the Hercules 
dSph, we adopted the specially fit 'colours', c\ and c-i, 
which represent a combination of our three filter photom- 
etry in B, V and r. These principal colours are designed 
to enhance the Hercules-to-field contra st in the CMP (see 
jSt, a nd are similar to those used by lOdenkirchen et all 
(2001). By fitting a ridge- line locus in this plane using 
the stars at the centre of Hercules, we found the prin- 
cipal colours to be: c x = 0.944(5 - V) + 0.330(1^ - r), 
c 2 = -0.330(B -V) + 0.944(V - r). A diagram in the 
(C2 , V) plane is essentially the CMD seen edge-on, and 
the dispersion of Hercules member stars around ci = is 
defined by the photometric errors. For the remainder of 
this paper, our CMDs will be defined in the (c\, V) plane. 
The blue errorbars in Fig. Q] represent the combined un- 
certainty of the colour ci, and the dashed lines represent 
the magnitude at which the photometric completeness has 
fallen to 50%: B50 = 25.60 ±0.05 mag, V 5Q = 25.72 ±0.05 
mag, r 50 = 25.56 ± 0.04 mag. 

3. THE STRUCTURE OF HERCULES 



The CMD of all stellar sources in the Hercules core re- 
gion is shown in Fig. Q] the first data for this galaxy to 
reach well below the main sequence turnoff. From the B 
and V photometry we derived the best-fitting distance, age 
and metallicity for a single stellar population (Table [J) us- 
ing th e CMD-fitting techniques descr ibed bvlde Jong et al.l 
(120071) and originally developed by iDofphinl (|2002ft . We 
find Hercules to be dominated by old and metal-poor stars. 
The dashed line in Fig. Q] traces the isochrone for a stellar 
population with an [Fe/H] abundanc e of —2.26 and an age 
of 13 Gyr (|Girardi et alJl2002L l200l . A detailed analysis 
of the star formation history and metallicity evolution of 
Hercules will be presented in a later publication. 

Starting from the CMD of all detected sources in the 
entire field, we tried to isolate probable Hercules dSph 
stars by excising "field sources" , such as foreground stars 
and background galaxies, using a CMD se l ection technique 
originally introduced by iGrillmair et al.l (|l995l ): see also 
lOdenkirchen et aT1l2001ft : a comparison was made between 
the (ci , V) CMD of the central region of Hercules and that 
of the field region (a distant portion of the 23' x 23' field 
of view) to produce a 'signal-to-noise' map of the Hercules 
stellar population across the CMD. The CMDs for the core 
of Hercules and the field are shown in Fig. [1] where the 
associated spatial regions were chosen iteratively and are 
outlined in Fig. [2] To make a map of probable Hercules 
member stars (Fig. [5]), we then chose a limiting signal-to- 
noise in the CMD plane (s = 3.0, the red contour in Fig. 
[1]) which minimised the field population while maintaining 
a high number of Hercules stars. This corresponds to the 
main sequence turnoff and upper main sequence region of 
the Hercules CMD; only stars above the 50% completeness 
limit shown in Fig. [T] were selected. 

At such faint magnitudes, background galaxies form the 
majority of contaminating sources, consequently optimal 
star-galaxy separation is important. The FWHM of many 
background galaxies are expected to be significantly larger 
than that of a star; hence, the DAOPHOT 'sharpness' pa- 
rameter provided a secondary level of filtering. We used 
the sharpness values from the B-band image, for which the 
PSF shape was effectively constant across the field and 
the seeing was best (0.8"). We constructed a sharpness 
map across the LBT i?-band image using the bright stars 
(14 < B < 22) and removed all sources with a non-stellar 
sharpness. Together, the CMD selection and star-galaxy 
separation techniques removed 85% of all stars, most of 
them from the field population. 

The spatial distribution of the 'cleaned' Hercules dataset 
is shown in Fig. [2] (upper panel). A relatively bright star 
(g = 15.5; SDSS) near the centre of the Hercules dSph 
is saturated in our data. The data in this region were 
excluded from our analysis. We converted this map into 
a stellar surface density contour diagram (Fig. [2J lower 
panel) by convolving each star with a Gaussian of radius 
0.6' (independent of stellar magnitude), where the con- 
tours trace stellar densities of 1.5cr, 3c, . . . , 10.5<r above the 
field density, where a is the variance of the mean back- 
ground level within a circle of radius 0.6'. 

Both panels in Fig. [5] show a strikingly elongated stellar 
distribution of the Hercules dSph. We quantified the in- 
ternal structure of Hercules by fitting a series of radially 
increasing elliptical contours to the stellar density map dis- 
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played in Fig. fusing the IRAF routine ellipse. The best- 
fitting position angle, ellipticity and central coordinates 
were measured at semi-major radii of 0.25', 0.75', . . . , 9.75' 
for 200 random data subsets (50% of all stars), with asso- 
ciated uncertainties calculated using bootstrap-sampling. 
The results in Fig. [3] indicate that the ellipticity of Her- 
cules is relatively mild towards the central (e ~ 0.3), 
however this increases sharply to an ellipticity of e ~ 0.65 
in the outer regions. That is, the major-to-minor axis ra- 
tio of Hercules is approximately 3:1, flatter than any 
other known Milky Way dSph except the tidally disrupted 
nearby Sagittarius dSph. 

We also constructed an azimuthally averaged stellar 
density radial profile of the dSph in a series of concen- 
tric annuli (separated by 0.5' with a fixed ellipticity of 
e = 0.65 and position angle of —73°). The background 
level was calculated as the density of stars in the upper 
(AS > 8') and lower (AS < -5') regions of the LBT field. 
This was subtracted from the profile datapoints, an d the 
result is shown in Fig. [3] We derived a best-fitting iKinel 
(1962) profile for this system using bootstrap-resampling, 
yielding a half-light (major axis) radius of 4.37 ± 0.29 ar- 
cmin (168 ± 11 pc), a core radius of 4.74 ± 0.57 arcmin 
(182 ± 22 pc), and a tidal radius of 25.9 ± 11.1 arcmin 
(994 ± 426 pc). Note that the formal tidal radius of the 
Hercules dSph is beyond the spatial range of our dataset, 
and should be treated with caution. 

4. DISCUSSION 

In their discovery paper, iBelokurov et alj (|2007t ) exam- 
ined the Hercules structure based on photometry complete 
to i ~ 22.5 (the depth of the sub-giant branch) and found 
that the system has an extended elongated morphology. In 
our analysis, we have probed ^3 magnitudes deeper, allow- 
ing a structural map of much greater signal-to-noise. Our 
data clearly confirm that this object is highly elongated 
(3 : 1) and shows some tentative indication of clumpy sub- 
structureQ (see Fig. [2]). 

It is possible that the elongated structure of Hercules is 
intrinsic to this system. If so, it would be the most flat- 
tened of all the Milky Way dSphs (except the Sagittarius 
dSph). All other known stellar systems with 3 : 1 flatten- 
ing or more are predominately supported by rotation; yet, 
no stellar system with a characteristic velocity as low as 
Hercules (a ~ 5 km s _1 ; ISimon fc Gehal [20071 is known 
to be rapidly rotating. If in equilibrium, Hercules would 
have unique structural or kinematic properties. 

Structurally, the Hercules dSph most closely resembles 
the Ursa Minor dSph (e ~ 56), a Milky Way satellite that 
is presumably dis rupting (|Martmez- Delgado et all l2001t 
iPalma et ai]|2003[ ). The only other highly elliptical Milky 
Way dSph is the Sagittarius system, a relatively bright 
dSph that is being disrupted by Galactic tidal forces. Tidal 
forces could therefore be a plausible explanation for the en- 
longated Hercules structure. However, the large distance 
of Hercules is problematic for this interpretation. Tidal 
forces decrease as i? -3 , and the presence of tidal distor- 
tion would imply that Hercules is on a highly elliptical 

11 Our measurement for the ellipticity of the inner regions is likely to be 
function described above. Also, the ellipticity is measured in the plane 
sight. 

12 Note that the visual impression of sub-structure in Fig. [2] is affected 
hashed region. 



orbit around the Milky Way, to allow the system to expe- 
rience strong tidal forces during its pericentric passage. 

A simple condition for tidal disruption of Hercules at 
its orbital pericentre, i? per i, is ^Here/mere ~ ^Mw/^peri, 
where CHcrc and ctmw are the velocity dispersions of the 
Hercules dSph and the Milky Way respectively, rnorc is the 
limiting radius of Hercules. For rncrc ~ 250 pc (that is, the 
geometri c mean of the major a nd minor axes), fTH crc ~ 5.1 
km s _1 ([Simon fc Gehal [20071 ). and ctmw ~ 150 km s _1 , 
we find that only i? pCT i ~ 8 kpc would lead to disruption. 
Given that Hercules is now located at a distance of 132 kpc 
(and may not yet be at apogalacticon) this would imply 
an orbit for this system of extreme ellipticity, e > 0.90. A 
detailed analysis of the feasibility and plausibility of this 
scenario is warranted, but beyond the scope of this letter. 
At any rate, regardless of whether Hercules' highly flat- 
tened structure is due to tides, which appears on balance 
most plausible, or intrinsic to an equilibrium system, its 
properties are exceptional. 
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a lower limit due to the 'blurring' effect of the Gaussian convolution 
of the sky, and Hercules may be further elongated along the line of 
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Fig. 1. — Colour-magnitude diagrams for the central region of 
Hercules (left panel) and a field region of equivalent area on the sky 
(right panel). The photometry has been corrected for extinction and 
all sources outside the sharpness limit have been removed. Our prin- 
cipal colour, ci, has been designed to enhance the Hercules-to-field 
contrast in the colour-magnitude plane. The red das hed line is the 
isochrone o f a metal-poor stellar population at 13 Gyr (IGirardi et al.l 
2002, 2004). The blue errorbars in the left panel represent the un- 
certainty in the colour ci returned by the artificial star tests, and 
the dashed lines trace the V magnitude at which the photometric 
completeness has fallen to 50%. The red contour line traces the 
CMD-selection region for the Hercules main sequence population. 
Note the increased background galaxy contamination in the right 
panel at magnitudes fainter than V = 24. 



Fig. 2. — Upper panel: Spatial distribution of the Hercules CMD- 
selected objects, where the dashed line marks the limit of our LBT 
data. The blue shaded region represents the saturated star aligned 
towards the centre of the dSph. The solid red ellipse marks the core 
radius from our best-fitting King model (see Fig. |3j and the dashed 
ellipse outlines the inner border of the field region. Both the core 
and field populations were used to derive the CMD-selection limit 
shown in Fig. [T] These ellipses have an ellipticity of 0.65 and semi- 
major axis radii of 4.42' and 11' respectively. Lower panel: Contour 
diagram of the CMD-selected sources. Each star has been convolved 
with a Gaussian of width 0.6' arcmin. The contours correspond to 
stellar densities of 1.5c, 3<r, . . . , 10.5<r above the background, where 
a is the uncertainty in the background stellar density from Poisson 
statistics. 
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Fig. 3. — Radial profile of Hercules, determined from the CMD- 
selected dataset. The stellar density was evaluated within ellipses at 
every 1' major axis radius using an ellipticity of 0.65 and a position 
angle of —73°. The dashed line represents the best-fitting King pro- 
file, and the parameters with their associated bootstrap uncertain- 
ties are listed. The background level of 15.43 ± 0.28 stars/arcmin 2 
(Poisson statistical uncertainty) has been subtracted from all data- 
points. Inset panel: Ellipticity as a function of radius. Uncertainties 
were determined using bootstrap resampling. The central coordi- 
nates of this system demonstrate a mild dependence on radius, with 
a variation of ~0.3' in RA and ~0.2' in Dec over the radial range. 
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Table 1 

Properties of the Hercules dSph 



Parameter 


Value 


R.A. (J2000) 


16:31:02.0 


Decl. (J2000) 


12:47:29.6 


E(B - V) (mag) 


0.055 ±0.005 a 


(to - M) (mag) 


20.6 ±0.2 


Distance (kpc) 


132 ± 12 


[Fe/H] 


-2.1 ±0.2 


Age (Gyr) 


13 ±3 


King rh 


4.37' ±0.29' (168 ± 11 pc) 


King r c 


4.74' ±0.57' (182 ±22 pc) 


King r t 


25.9' ±11.1' (994 ±426 pc) 


c = log (r t /r c ) 


0.74 ±0.25 



"From Schlegel et al. (1998). The uncertainty quoted here represents the variation in reddening over the LBT field. 



